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Abstract. The absolute V light curves of 5 SNe la, se- 
lected to represent the known range of absolute luminosi- 
ties at maximum for this class of objects, are presented. 
Comparison of the long term luminosity evolution shows 



SNe la has shown the existence of a range in the proper- 
ties of this type of SNe. In particular, a correlation has 
been found between the absolute magnitude at maximum 



and the shape of the light curve (Hamuy et al., 1995; Riess 



that the differences seen at maximum persist, and actually ct al., 1995). If this relationship can be well calibrated and 



increase with time, reinforcing the notion that intrinsic 
differences do exist among SNe la. Since such differences 
are not accounted for in the standard progenitor scenario, 
it becomes important to derive constraints for the mod- 
els directly from the observations. In order to investigate 
the influence of the two most important parameters, that 
is the masses of the synthesized radioactive material and 
of the ejecta, a simple MC light curve model was used 
to simulate the luminosity evolution from the explosion 
to very late epochs 1000 days). It was found that the 
observations require a range of a factor 10 in the masses 
of the radioactive material synthesized in the explosion 
(A^Ni — 0.1 — I.IM0,) and a factor 2 in the total mass of 
the ejecta (Mej = 0.7 — I.4M0). Differences of a factor 2 
in Mni seem to be present even among 'normal' SNe la. 

Some evidence was also found that the deposition of 
the positrons from Co decay varies from object to ob- 
ject, and with time. In particular, the latest HST obser- 
vations of SN 1992A seem to imply complete trapping of 
the positrons. 



Key words: supernovae: general - supernovae: individual: 
SN 1991bg, SN 1991T, SN 1992A, SN 1993L, SN 1994D 



1. Introduction 

For many years type la supernovae (SNe la) were consid- 
ered very reliable standard candles. However, more accu- 
rate CCD photometry for a rapidly growing number of 



other parameters are found to play a negligible role, the 
confidence in the use of SNe la as distance indicators may 
be restored. 

A more subtle problem concerns the progenitors of 
SNe la. The standard scenario for SNe la is that a degen- 
erate white dwarf in a binary system accretes matter from 
a close companion, either red giant or white dwarf (WD), 
until it reaches the Chandrasekhar mass (Mph = 1.4Mq ^ 



and explodes, leaving no r emnant ( [Nomoto et al., 1984 ; 
Woosley and Weaver, 198^ ) . By tuning the relative masses 
and the separation of the stars in the binary system it is 
possible to obtain a very slow evolution and thus to ex- 
plain the occurrence of SNe la among population II stars 
and hence in elliptical galaxies. If all SNe la events are 
triggered by the WD's mass reaching Mch, independently 
of the initial conditions, it is natural to expect some ho- 
mogeneity in the outcome of the explosions. One should 
nevertheless keep in mind the fact that type la super- 
novae occur in all types of galaxies. This theoretical view 
reinforced the empirical use of SNe la as standard can- 
dles. Lately, however, with the advent of larger samples, 
the observed diversity among SNe la appears to challenge 
that paradigm. 

In particular, in the following we will show how the 
simple comparison of the late light curves of different SNe la 
requires the existence of a range in both the mass of ra- 
dioactive material synthesized in the SN la explosions and 
the mass of the ejecta. Such differences can be seen not 
only in the well-known cases of peculiar SNe such as the 
bright SN 1991T and the faint SN 1991bg, but also within 
the group that has hitherto preserved the definition of 
'normal' SN la. 



* Based on ESO observations collected at ESO-La Silla 
(Chile) 
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2. Late light curves 

We have compiled the absolute V light curves of five SNe la, 
namely SN 1991T, SN 1991bg, SN 1992A, SN 1993L and 
SN 1994D, combining unpubHshed observations with data 
from the literature. 

In principle, bolometric light curves should be used for 
comparison with the models. Unfortunately, such data are 
available only for very few SNe, even in the best case they 
are limited to the UVOIR spectral region and, with the 
exception of SN 1987A, they do not cover late epochs. 
However, both observations and spectrum synthesis indi- 
cate that in SNe la most of the deposited radioactive en- 
ergy is re-radiated in the optical region and that, at phases 
later than lOOd, the bolometric correction is constant and 
small, mboi-mv ^ 0-1 ~ 0-2 mag ( ^untzeff, 1996| ; |Wheeler| 
and Hiiflich, 1997). Therefore, in the following we simply 



assume that the bolometric correction for the V band is 
zero for all SNe la. Despite this crude approximation, this 
is not the major source of error, the uncertainties in the 
parent galaxies distance moduli and extinctions being in 
fact much larger. 

The observations were retrieved from the archive of 



the ESO supernova monitoring programme ( Turatto et al 
1990a|]| and are presented here for the first time, except 
for SN 1991bg, whose complete light curve has already 
been published in Turatto et al. ( 1996 ). These data have 
been supplemented, for the early phases, with published 
photometr y from Phillips et al. ( |1992| ) for SN 1991T, Sun- 
zeff et al. ( |1996D for S N 199 2A, Filippenko et al. ( |1992D 
and Leibundgut et al. ( |1993| ) for SN 1991bg and Patat ct 
al. ( |1996| ) for SN1994D. 

Additionally, we included the photometry of SN 1992A 
obtained with WFPC2 on HST on Aug 2, 1994, 926 days 
after maximum. These observations, obtained in the frame- 
work of the SINS program (Kirshner et al. 1993) , consist 
of two sequences of four exposures through the F555W 
and the F439W filters, which are similar to the V and 
B bands, respectively. Exposure times for the individual 
frames were 900 sec for F555W and 1200 sec for F439W. 
The individual frames, calibrated in the standard pipeline, 
have been retrieved from the HST archive, properly aligned 
and combined to eliminate cosmic rays. In the combined 
F555W image we found a stellar object whose offset from 
the field stars agrees to within 0.2 arcsec in both coordi- 
nates with the offset of SN 1992A as measured on the last 
ground-based observations (note that 1.2 arcsec west of 
the SN is a 23 mag background galaxy which appears as 
a stellar object on ground-based images). The SN magni- 
tude, as measured by means of aperture photometry and 
calibrated relative to a sequence of local stars whose mag- 
nitudes were determined from the ground-based observa- 
tions, is V — 25.9 ± 0.3. The SN is not so evident in the 
combined F439W image which, however, is consistent with 
a {B — V) colour close to 0. 



The main data for the five SNe la are listed in Tab. ^ 
while the absolute V light curves are shown in Fig. ||. Ab- 
solute magnitudes were computed using the distance mod- 
uh given in the TuUy ( |1988| ) catalog, except for SN 1993L, 
whose parent galaxy is not listed there. For this SN we 
use the distance given in the Lcda extragalactic database^ 
Uncertainty in the absolute magnitudes arises mostly from 
the uncertainty concerning the distance scale, (both cat- 
alogs we are using adopt a Hubble constant of 75 km s~^ 
Mpc~^). Although the errors in the absolute magnitudes 
may be quite large (of the order of ±1 mag), the rela- 
tive distances should be more reliable (typical errors are 
±0.2 — 0.3 mag) and so the differences in magnitude be- 
tween different SNe are expected to be real. 

Magnitudes have been corrected for the total absorp- 
tion Av as estimated in the references indicated in Tab. l[ 

The reference epochs for the light curves were chosen 
using the times of explosion estimated from spectrum syn- 
thesis and light curve modelling of the early photospheric 
phase. These range between 12 and 20 days prior to B 
maximum. The uncertainty in the determination of the 
time of maximum is negligible for our discussion, which is 
mostly concerned with the very late epochs. 

In Tab. |l| we also list the quantities Amfs, the dif- 
ference in B magnitude from maximum to 15d, which is 
widely used to characterize the early light curve (col. 8), 
and Am^QQ , the difference in V magnitude from maximum 
to 300d (col. 9). When observations at this epoch were 
not available, the value of Amg^Q was derived by interpo- 
lation between the closest adjacent observations. Finally, 
we report the maximum expansion velocities of the Fe- 
nebula (col. 10). These are derived from models of the 
emission lines in spectra at epochs around 300d, except 
for SN 1991bg, whose spectral evolution was fast and for 
which a spectrum at 221d was used (Mazzali et al., 1997). 
The velocities tabulated correspond to the outer veloc- 
ity of the Fe-nebula for which synthetic nebular spectra 
give a best fit to the late-time spectra of the SNe at hand 
(Mazzali et al., in preparation). 

The SNe have been selected to represent the known 
range of luminosities of SNe la at maximum, going from 
the faint SN 1991bg to the 'average' SNe 1992A, 1993L 
and 1994D and the bright SN 1991T. The fuU range is 
about 2.3 mag at maximum. In particular, it appears that 
SN 1994D is consistently and significantly brighter than 
SNe 19 92A and 1993L. This was also implied by Hamuy 
et al. ( 1995 ) on the basis of the different decline rates. 
Fig.l shows that the differences in absolute magnitude 
persist to the very late phases, and actually increase with 
time, mostly reflecting the range in the early decline rates. 
Differences are about 3.5 mag already 2 weeks after maxi- 



^ The Lyon-Meudon Extragalactic Database (LEDA) is sup- 
plied by the LEDA team at the CRAL-Observatoire de Lyon 
(France) . 
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Table 1. SNe la data 



SN 


Gal 






M ~m 


Av 


Mv 






Neb.vel. 

[km s"^] 


1991bg 


NGC4374 


13.96 


0.74 


31.13 


0.15" 


-17.27 


1.95" 


8.4 


2500 


1993L 


IC5270 


13.2* 


0.2 


31.52 


0.23* 


-18.5# 


1.5# 


7.4 


8000 


1992A 


NGC1380 


12.54 


0.00 


31.14 


0.00'' 


-18.60 


1.47"* 


7.3 


8500 


1994D 


NGC4526 


11.90 


-0.08 


31.13 


0.08^= 


-19.31 


1.26= 


7.3 


9500 


1991T 


NGC4527 


11.50 


0.13 


30.65 


0.40-^ 


-19.55 


0.94s 


6.7 


10500 



* - measured at the epoch of the B maximum 
$ - only galactic absorption included 

# - the SN was discovere d ten days after maximum. Its m agnitude at maximum was ex trapolated by comp arison with SN 1992A. 
Refs.: a - Turatto et al. ( [l996|); b - Kirshner e t al. ( |l99j ); c - Ho & Fihppenko ( |l995| ); d - Hamuy et al ( [l996| ); e - Patat et al. 
([1996|); / - Mazzali et al. (|1995|); g - Phillips (|l993|) 



mum, j.nd reach about 4 mag at 300d (Tab. l[), after which 
they apparently stop increasing. 

The exception is SN 1991T, for which the huiiinos- 
ity after 500d stops declining. This results from an echo 
formed as the light emitted by the supernova near max- 
imum light was reflected off dust in the CSM or ISM. 
(Schmidt et al., 1994; Danziger et al. unpublished spec- 
tra). Since there is no direct relation between the late- 
time luminosity of SN 1991T and the SN remnant itself, 
in the following we ignore the observations of SN 199 IT 
at phases later that 500d. 

^From Tab. |l| the correlation between the absolute 
magnitudes at maximum and the early decline rate {Amf^) 
is evident. In addition, it appears that at late phases 
brighter SNe have larger Fe-nebula expansion velocities 
(as already noticed by Danziger, 1994 and Turatto et al., 
1996 ) and slower luminosity decline rates. 

Another interesting feature shown in Fig. ^ is that, 
although the late light curves may appear quite linear 
on short time intervals, the decline rate actually slows 
down in the long run (> bOOd). Unfortunately, the ev- 



et al., 1990b), implying the presence of a delayed energy 



idence lis based on only one observation tor each ot the 
SNe igyibg, 1992A and, possibly, I993L. Since obscrva- 
tions at such late epochs are difficult, and the photometry 
may be contaminated by phenomena not directly related 
to the remnant itself (eg. echoes in the case of SN 1991T 
or contaminating unresolved stars) this suggestion needs 
further confirmation. Meanwhile, it is worth exploring the 
implications of the effect if real. 



3. Ra dioactive deca^r 



input. The linear tail of the light curves corresponds to 
an exponential luminosity decline, suggesting that the ra- 
dioactive decay of unstable heavy elements is a likely en- 
ergy source, even though the decline rates for SNe la do 
not match the decay rates of known radioactive species. 

Indeed, model calculations have shown that explosive 
silicon burning can produce unstable isotopes of iron group 
elements, in particular ^^Ni. ^^Ni decays into ^^Co, emit- 
ting 7-rays with an average energy per decay of 1.71 MeV 
( [Burrows and The, 1990[ ). Since the half-life of ^^Ni is rela- 
tively short (~ 6.1 days), this decay can only be important 
for the early light curve (the first 1-2 months). 

^^Co is also unstable and decays into '"^^Fe. The aver- 
age energy available per decay is 3.67 MeV. Although most 
of this energy is released in the form of 7-rays, a signifi- 
cant fraction of the ^^Co decays (19%) produces positrons. 
Positrons deposit their kinetic energy in the ejecta and 
then annihilate with electrons, producing two photons of 
energy E-y = irieC^ each. The positron kinetic energy ac- 
counts for a bout 3.5% of t he total ^^Co decay energy (Ar 



nctt, 1979| ; Axelrod, 1980 ). Because of its longer half- life 
(77.7 days), ^^Co decay can power the light curve of SNe 
for at least 2-3 years. 

There is now ample evidence that the radioactive decay 
of ^^Ni and ^^Co provides most of the energy during the 
first 2-3 years, as was first suggested by Pankey (1962) and 
then, independently and more quantitatively, by Colgate 
and McKee (1969). The most convincing direct evidence 
has been the detection of 7-rays lines from ^^Co decay in 



the t y pe II SN 1987A (pVIatz et al., 1988| ; [Arnett et al , 
1989| ; Palmer et al., 1993| ) a nd the temporal beha viour 



Starting with the earliest attempts to model the SN light 
curves it has been clear that the fast expansion of the 
ejecta would result in a very rapid cooling and hence in a 
very rapid luminosity decline if only the thermal energy 
from the explosion had been available to power the SN. 

On the contrary, for almost all types of SNe, the lumi- 
nosity evolution is relatively slow. In particular, at epochs 
later than 150-200d the light curves decline almost lin- 
early with rates in the range 0.8-1.5 mag/lOOd (Turatto 



of the [ Coll] lines at 10.52 /i ( Danziger et al., 199C| ) and 
1.547/i ([Varani et al., 1990| ). In the case of many SNe la, 
modelling of the late-time spectra indicates the presence 
of ^^Co lines whose intensity evolution is in agreement 
with the prediction f rom the ^^Co decay ( [Axelrod, 1980 ; 



Kuchncr et al., 1994) 



In the early phases of the SN evolution the density of 
the ejecta is still high enough that the 7-rays from radioac- 
tive decay are trapped in the ejecta and completely ther- 
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200 400 600 800 1000 
days from explosion 



Fig. 1. The absolute V light curves of 5 SN la are compared with 4 models with different ejecta and Ni masses. The 
reference epochs is the time of explosion estimated from spectrum synthesis and light curve modelling of the early 
photospheric phase. From top to bottom the models are characterized by M^j = 1.4, Afjvi — 1.1 (continuous line); 1.4, 
0.8 (dotted); 1.0, 0.4 (short dashed) and 0.7, 0.1 (long dashed). For all models k+ = 7 is adopted. In the inset, the 
epoch near maximum is shown enlarged. 



malized through Compton scattering with the free elec- 
trons. With the expansion the density decreases and the 
mean free path of the 7-rays increases. In the case of SNe II 
(and of at least a fraction of the SNe Ib/c) the mass of the 
ejecta is so large (> 3Mq) that even if their density de- 
creases they remain optically thick to 7-rays for 2-3 years 
(Woosley, 1988). On the other hand, the ejecta become 



transparent to optical radiation a few months after max- 
imum, and the observed decline of the optical luminos- 
ity reflects directly the decline of the radioactive energy 
supply. Indeed, the exponential tails of most SNe II light 
curves show the same e- folding time as that of the ^^Co 
energy release. 

In the case of SNe la, the ejecta are less massive, while 
the Ni mass is larger (AfNi/Afcj ~ 0.5 in SNe la vs. 0.01 
in SNe II). Also, there is evidence that in some cases Ni 
can be found in the outer layers. This allows an increasing 



fraction of the 7-rays to escape thermalization as time 
goes on, so that the light curves of SNe la decline more 
rapidly than the ^^Co energy release. 



4. A simple light curve model 

The observed differences in the light curves of SNe la may 
be attributed to many different factors. Accurate fitting of 
the light curve of a particular SN la may indeed require an 
ad-hoc explosion model and fine tuning of several param- 
eters. However, this approach makes it difficult to disen- 
tangle the relative importance of the different parameters. 
Since we wanted to isolate the effects of the two arguably 
most interesting physical parameters, Mni and Moj, we 
developed a simple Monte Carlo code for the calculation 
of the bolometric light curve and did not use any partic- 
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ular explosion model. In particular, we made no effort to 
fit the shape of the light curve near maximum in detail. 

Starting from a given density stratification and a dis- 
tribution of Ni within the ejecta, the code computes the 
deposition of the 7-rays as a function of time with an MC 
scheme. The 7-ray energy sources arc given by Sutherland 
& Wheeler (1984), while the opacity to 7-rays is assumed 
to be gray and to have a constant value = 0.027cm^ 



g ^ , in agreement with recent calculations ( Swartz et al 
1995|)~|The energy released by the decay of ''''Co in the 
form of positrons (19% of the total) can only be trans- 
formed into 7-rays if the positrons annihilate. 

The calculations show that with the expansion the 
ejecta become rapidly transparent to 7-rays. In particular, 
200 days after the explosion the fraction of the 7-ray en- 
ergy which is deposited in a typical model is only ^ 0.5%. 
At this time the kinetic energy of the positrons may be- 
come the main contributor to the light curve, if a major 
fraction of it is deposited. 

The early assumption was that positrons deposit all of 



their kinetic energy and annihilate almost on the spot ( Ax- 
elrod, ,980). Subsequently, this idea has been questioned. 



For instance, Chan & Lingenfelter ( 1993 ) suggested that 
a significant fraction of the positrons manage to escape, 
possibly because of the presence of radially combed mag- 
netic fields. On the other hand, Swartz et al ( |1995[ ) argued 
that positrons cannot escape at all, otherwise a significant 
fraction of the radioactive decay energy would be lost. 

Recently, it has been argued that to fit the late light 
curve of SNe la requires that the ejecta become progres- 
sively transparent also to positrons ( [Colgate et al., 1997 ) , 
although with a longer time scale than for 7-rays (a typ- 
ically adopted value is = 7 ciyP g^^)- Describing the 
positron deposition by means of a positron opacity ap- 
pears to be a natural approach, and we have adopted it 
in our light curve code. This gives us the freedom to ma- 
nipulate the positron deposition as well as the transport 
of optical photons. 

Our 7-ray source function for '^^Co therefore takes the 
form 

Sco^0Ms + 0.19sD^+ (1) 

where s = 6.78 10^[exp(— t/rco)— exp(— t/rAr;)] ergg^-^ s^"^ 
is the rate of energy production from the ^®Co decay (ex- 
cluding the KE of the positrons), tco and r^i are the 
mean lifetimes of ^^Ni and ^^Co, respectively, and Dg+ 
is the positron deposition function. The positron KE ac- 
counts for an additional 7-ray source, which is given by 

S^+=0m6sD^+. (2) 

When the 7-rays and the positrons deposit their en- 
ergy, it is assumed that after thermalization this energy 
emerges as optical photons, so the contribution to PdV 
work is negligible. These optical photons can only con- 
tribute to the bolometric luminosity when they manage 
to escape from the ejecta. Given the high ejecta density, 
especially early-on, escape is not instantaneous. This is 



the reason why SNe la reach their LboI peak only 2-3 
weeks after the explosion. The random walk of the optical 
photons in the expanding ejecta is also followed with an 
MC scheme. A gray optical opacity Kopt is assumed where 
only a pure scattering opacity obtains. The ejecta are as- 
sumed to be purely scattering. Time delay is accounted for 
by computing the time elapsed between successive scatter- 
ings, and updating the density before a new free flight is 
started. So finally for each optical energy packet we have a 
'production time', corresponding to the time of emission, 
i.e. the time of deposition of the 7-ray packet from which 
the optical packet originates, and an 'emission time', the 
time when the packet eventually escapes from the ejecta. 
The difference between these two times becomes smaller 
and smaller as time goes on. The code allows us to repro- 
duce not only the declining part of the light curve but also, 
at least roughly, the rising branch and the near-maximum 
phase. For 'average' SNe la we find that Kopt = 0.15 cm^ 
g~i gives a reasonable fit to the rising branch and to the 
maximum of the light curve, whereas the late light curve 
is quite insensitive to the value of this parameter. 

With these assumptions, the most important param- 
eters for the light curve modelling are the mass of syn- 
thesized ^^Ni (MNi), which determines the absolute lu- 
minosity, and the total mass of the ejecta (Moj), which 
determines the optical depth for the radiation from the 
radioactive decay. Hopt affects mostly the time of maxi- 
mum and its brightness, while infiuences the late-time 
behaviour. For larger Kopt the maximum occurs later, and 
it is fainter and broader. For smaller Kp+ the decay at late 
times is faster. 

In the past, the consensus was that M^i and M^j were 
exactly the same in all SNe la events, making them ideal 
standard candles. Estimates for A/ni were in the range 0.5 
to 0.7 Mq, while M^j = Mch- However, there is increasing 
evidence that both Mni and Mej can be very different in 
different objects, although most SNe la appear to cluster 
around the historical values. 

Since the SNe we have listed in Table 1 span a range 
of absolute magnitudes at maximum and of decline rates, 
it is quite likely that they cover a range of AfNi, and pos- 
sibly also of Mcj. To investigate this question further, we 
used our MC light curve code. As we discussed above, Mni 
and Afcj can be easily changed as input to the code. We 
took the approach of using a simple code and of not us- 
ing any particular explosion model for any particular ob- 
ject, because literally hundreds of models, involving differ- 
ent Afcj and explosion mechanisms, and producing differ- 
ent amounts of ^^Ni, have now been produced by various 
groups, making it difficult to distinguish between them on 
the basis of their ability to fit observed light curves. 

We adopted a W7 density structure, and rescaled it 
according to the epoch and to the values of A/ni and Mgy 
We assumed that the Ni resides at the centre of the ejecta, 
in rough agreement with models of SN la explosions. The 
only exception is the model for SN 1991T. In this case. 
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analysis of the early-time spectra suggests that there is 
an outer Ni shell (Mazzali et al. 1995). 

Homologous expansion is assumed for rescaling to the 
appropriate epoch. Since different values of A/n; and M^j 
should lead to different kinetic energies of the ejecta, the 
p{v) dependence is also rescaled. If we define the ratios 
X^i — AfNi/0.60Mo and X^j = Mcj/Mch, we can rescale 
the kinetic energy per unit mass according to KE oc 

In doing this we assume that all the KE is produced by 
burning to ^^Ni, neglecting the contribution of incomplete 
burning to intermediate mass elements (IME) such as Si. 
This is a reasonable approximation as long as A//j\/ e /A/ni 
is similar for all values of Afwi, which is of course not neces- 
sarily true. In particular, the energy will depend not just 
on A/ni, but also on a third parameter, the total mass 
of the newly synthesized elements. Indeed, in some cases, 
like SN 1991bg, the mass of these elements (e.g. Si-Ca, 
^'^Fe and ^^Ni) is comparable to, and probably greater 
than AfNi. Thus we expect our models to be an increas- 
ingly poor representation of the real properties of the SN 
the further A/ni departs from the reference W7 value of 
O.6M0. 

On the other hand, most SN la explosion models, in- 
cluding models for sub-Chandrasekhar explosions (Woosley 
& Weaver 1994), produce p{v) functions whose shape is 
similar to that of W7, so in this respect our approximation 
should be reasonable. The velocity of a given shell in the 
ejecta is then rescaled according to u/t)^7 — {XNi/Xej}^^^ , 
while the density of that shell is given by p/pw7 — 

y5/2 t^-3/2 

Note that models with a ratio AfNi/Aicj — 0.43 have 
the same KE per unit mass as W7 (for which A^Ni/A^cj = 

0. 60/1.38), and their density is simply the W7 density 
structure times X^j , while models for which X^j'^ — X^j^ 
have the same density as W7, but with each shell shifted 
to a velocity v = vwii^jf-^Y^"^ ■ Since in the photosphcric 
epoch the apparent photosphere forms at an almost con- 
stant value of the density, this may be a useful indepen- 
dent tool to investigate the values of A^Ni and Afoj. This 
approach was adopted for all sub-Chandrasekhar models, 
but not for the SN 1991T ones, since the observed outer 
distribution of ^^Ni clearly must have a different effect on 
the kinetic energy. Thus, the KE for the SN 1991T model 
is probably underestimated. 

Apart from the somewhat arbitrary rescaling of the 
velocity field, several other uncertainties must be kept in 
mind when comparing the model bolomctric light curve 
with the observations: 

1. Though we argued the V and bolometric light curves 
are not very different (cf. Sect. 2), they are not the 
same. A proper comparison would require a full NLTE 
light curve calculation, with the appropriate SN model. 
Such calculations do not exist as yet. Actually, even the 
bolometric light curve would require a more complex 



calculation than in our simple code. In particular, the 
opacities may change as a function of depth in the 
ejecta, and with time. 

The SN may be clumpy, thus changing the deposition 
rate. This would also influence the nebular spectrum 
(Mazzali et al., in preparation). 

The positon deposition is not well known. We discuss 
this point in the next section. 

The steady-state assumption that the luminosity equals 
the instantaneous energy input breaks down after about 
600 days, when the deposition time for the primary 
electrons produced by the 7-rays becomes long with 
respect to the dynamical and ionization time scales of 



the nebula (Axelrod, 1980). The ionization state may 
then be far from a steady state solution. 
5. The distance and extinction to the various SNe are un- 
certain. By using a single source for the distances, we 
should at least minimize the relative errors, but the ab- 
solute numbers could change. For instance a Cepheid- 
based set of distances would require much bigger Ni 
masses. 

5. Results 

We selected the model parameters in order to fit the mag- 
nitude at maximum and the long-term decline rate of the 
individual SNe. Although less emphasis was given to the 
detailed fit of the light curves, we comment on the major 
discrepancies and the possible reasons. 

In a first set of models we kept ~ 7 constant, and 
tried to obtain the best fit by varying A/ni and Afcj. In 
general, increasing A/ni causes a shift of the light curve to 
brighter absolute magnitudes, while increasing Afcj delays 
the time of maximum and makes the decline of the light 
curve less steep. 

The model fitting to the V light curves of the five 
SNe la is shown in Fig. |l|. The models used for the fits 
are as follows. 

For SN 1991bg, A/ej = O.7M0, Mm = O.IMq. This 
is in agreement with findings from modelling of both the 
photosphcric and the nebular-epoch spectra (Mazzali et 
al. 1997), making SN 1991bg the best-studied case of a 
probable sub-Chandrasekhar mass SN la. 

The V light curves of SNe 1992A and 1993L are prac- 
tically indistinguishable once they have been corrected 
for extinction and relative distance. With the assumed 
distances, they reach a maximum of about —18.6 mag. 



The curves are fitted by a model with A/o 



l.OM 



AfNi = 0.4Mq, which suggests they may also be sub- 
Chandrasekhar events. The model has the same AfNi/Afoj 
ratio as W7. 

SN 1994D appears to be brighter than SNe 1992A 
and 1993L. The peak magnitude (—19.4) is well fitted 
by a model with Afoj = 1.4M0, M-^i = O.8M0, imply- 
ing that the Ni production may be higher than in 'clas- 
sical' SNe la. That SN1992A is fainter than average has 
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also been demonstrated recently by Delia Valle et al. (in 
preparation) who determined the distance to its parent 
galaxy using the method of the globular cluster luminos- 
ity function. 

As can be seen from the insert of Fig. 0, a common fea- 
ture of the models is that the rise to maximum is steeper 
and the early decline slower than observed. This could be 
improved by adopting different Ni distributions. In gen- 
eral, placing Ni further out than in the centre of the ejecta 
makes the light curve faster. We have computed a model 
using the W7 Ni distribution, which has a 'hole' in the 
central, highest density regions. This leads to faster es- 
cape of the optical photons after maximum, and produces 
a faster-declining light curve, which compares well with 
the observed ones. Still, W7 seems to be too faint to re- 
produce the light curve of SN 1994D and too bright for 
that of SN 1992A. 

It should be noted here that the V light curves of 
SNe 1994D, 1992A and 1993L are very similar in shape. 
Thus, if the relative distances and extinctions were appro- 
priately adjusted, all three SNe could be explained with 
a single explosion model. On the other hand, significant 
observational evidence also exists in favour of there being 
a difference between SNe 1994D and 1992A: SN 1994D is 
about 0.1 mag bluer at maximum, has a lower Sill line 
velocity near maximum (Patat et al. 1996, Fig. 10a, where 
SNe 1994D, 1989B and 1990N appear to form one group 
and SNe 1992A and 1981B another), and broader nebular 
lines, with the SNe just mentioned falling again into two 
different groupings (Mazzali et al, in prep.). Thus, intrin- 
sic differences between SNe generally labelled as 'classical' 
la are only beginning to receive the attention they deserve. 

The last object we tried to fit is SN 1991T. For this 
SN, analyses of the early- (Mazzali et al. 1995) and late- 
time (Spyromilio et al. 1992) spectra suggest that the Ni 
mass is about IMq, and that a significant fraction of it 
is located in the high velocity outer part of the ejecta. If 
we assume Moj = Mch, a good fit can be obtained for 
Mni = I.IM0, of which O.6M0 is in the centre and O.SMq 
is in the outer layers, which confirms previous findings. 

Thus, it appears that a range of almost one order of 
magnitude in Mni is required to fit all the objects, and of 
at least a factor 2 in Afcj- The range in Mni reflects the 
range in absolute magnitudes at maximum rather closely. 

The models with Kg+ = 7 fit the observations up to 
400-500 days well for all the objects except 1991bg. In the 
case of SN 1991bg, starting 100-150 days after explosion 
the model is brighter than the observations. To reconcile 
the model with the observations, one could further de- 
crease A/ej, but this would lead to a very early maximum 
and would also cause problems for the interpretation of 
the spectra near maximum (Mazzali et al., 1997). Actu- 
ally, a good fit to the late light curve of SN 1991bg can be 
obtained assuming that the opacity for positrons is much 
smaller than Kf,+ = 7. This is shown in the bottom panel 
of Fig. ||, where the model is calculated for the extreme 
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Fig. 2. Comparison of the V light curves of different 
SNe la with models for different values of the positron 
opacity, k+. The models for the individual SNe are char- 
acterized by different values of the radioactive Ni and 
ejecta masses. From top to bottom these are: 1991T: 



1.4Mq, Af, 



I.IMq; 1994D: M, 



lAMr. 



Mm = 0.8Mq; 1992A and 1993L: 
0.4Mq; 1991bg: A/oj = O.7M0, Mm = O.IM^ 



ej — -L.1J-V-L0, 

I.OMq, A^Ni = 
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case of complete transparency of the ejecta to positrons 

{Ke+ = 0). 

The opposite trend may be indicated by the SNe la 
1993L and 1992A. At phases later than 400 days the ob- 
served light curves appear to decline at a rate slower than 
the prediction of the model with k,^+ = 7. For SN 1994D, 
no observations are available at these very late epochs. 

Finally, the light curve of SN 1991T is compatible with 
Ke+ > 7. This may be real, but it may also indicate that 
the fraction of Ni on the outside is less than what we 
have assumed. Note, however, that if we distribute all the 
I.IM© of Ni in the centre and compute a rescaled model 
with Kg+ — 7, this has a larger KE than W7, and there- 
fore a lower density, so it actually declines faster than the 
Ke+ = 7 model shown in Fig. 2. Another possibility is that 
Mcj > Mch- We computed a model with M^] = l.GM©, 
-^^Ni — 1.2Mq, which produces a light curve with a broad 
maximum and a slow decline. This model fits the obser- 
vations reasonably well at all epochs, including the late 
phases. In the phase 60-150 days the model light curve 
is brighter than the observed one, but this is a feature 
common to all models shown in this paper. Even in this 
case, however, Kg+ > 7 cannot be ruled out at 300-400 
days. Clearly, SN 1991T deserves a much more detailed 
investigation than has been attempted here. 

In principle we cannot exclude the possibility that, in 
addition to the radioactive decay, something else may con- 
tribute to the luminosity. In particular, given that spec- 
tra are not available at these very late epochs, we cannot 
entirely rule out contributions from echoes by circumstel- 
lar dust, as has been observed for SN 1991T. However, 
the luminosity dechne of SN 1992A from 400 to 926 days 
(0.9 ± 0.1 mag/lOOd) is very close to the ^^Co decay rate 
(0.98 mag/lOOd,) and other energy sources do not seem 
to be required. As shown in Figj^, this is the decline rate 
expected in the case of essentially complete trapping of 
the positrons (Ke+ — 300) and complete transparency to 
the 7-rays. 

6. Conclusions 

The observed V light curves of different SNe la have been 
compared with models of the bolometric light curve ob- 
tained varying the values of Mm and A/oj . The distribution 
of Ni was assumed to be central with the exception of the 
models for the bright SN 1991T, where an outer Ni shell 
was also used. The KE of the ejecta was rescaled accord- 
ing to Mn! and Afoj in all cases except SN 1991T. Both 
Mni and Afoj infiuence the light curve, but while Mm in- 
fluences mostly the level of the curve, Afoj affects also the 
time of maximum and the decline rate, so that a simulta- 
neous fit to both the maximum and the late light curve, 
up to when the positrons begin to dominate, can only be 
obtained for a single set of Mm and Mcj values. Changes 
in the distribution of Ni and in the opacities between the 
various models would of course alter this result. 



The comparison of the (late) light curves of SNe la 
with models confirms the suggestions derived from the 
analysis of observed spectra near maximum of the exis- 
tence of a range in both the Ni and the ejecta mass. 

It appears that objects which have usually been re- 
garded as 'typical' SNe la may actually differ by more than 
0.5 mag at maximum. This difference increases with time, 
suggesting that both Mmi and Afoj may be significantly dif- 
ferent in these objects (in particular, SN 1994D appears 
to be brighter and to have a larger Afoj than SNe 1992A 
and 1993L). Not surprisingly there appears to be a corre- 
lation between Mmi and Afoj for the sample of 5 objects 
discussed here. 

In general a reasonable fit to the late light curves can 
be obtained by assuming that the ejecta become progres- 
sively transparent also to the positrons generated in the 
Co decay. However, there are indications that the positron 
opacity is not constant from object to object or with time: 
the late light curve of SN 1991bg can be fitted if the as- 
sumption of complete transparency of the ejecta to the 
positrons is made, while those of SNe 1992A and 1993L 
after 500 days seem to require complete trapping of the 
positrons. 
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